Synchrotron x-ray study of lattice vibrations in CdCr 2 04 
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Using inelastic x-ray scattering we have investigated lattice vibrations in a geometric frustrated 
system CdC^CU that upon cooling undergoes a spin-Peierls phase transition at Tn = 7.8 K from a 
cubic and paramagnetic to a tetragonal and Neel state. Phonon modes measured around Brillouin 
zone boundaries show energy shifts when the transition occurs. Our analysis shows that the shifting 
can be understood as the ordinary effects of the lowering of the crystal symmetry. 

PACS numbers: 63.20.-e, 78.70.Ck, 



A spinel AB2O4 system is an excellent model to study 
the physics of frustration. The octahedral B sites sur- 
rounded by oxygen ions form a three dimensional network 
of corner sharing tetrahedra, called a pyrochlore lattice. 
Since oxygen octahedra form an edge sharing network, if 
the B site is occupied by a magnetic ion with unpaired 
t2g electrons, the nearest neighbor interactions become 
dominant, which yields strong frustration. In spite of the 
theoretical predictions that the pyrochlore system should 
not order at any temperature, spinels usually undergo 
phase transitions into ordered states at nonzero temper- 
atures 1 -. This is because the spin degree of freedom can 
be coupled with other degrees of freedom, such as orbital 
and lattice, to lift the magnetic degeneracy 2 -. 

One example is a spin-Peiels phase transition that oc- 
curs in ACr204. The Cr-based spinels in the absence of 
the orbital degree of freedom, AC^CU, remain param- 
agnetic to temperatures far below the Curie- Weiss tem- 
perature, -390 K and -88 K for A = Zn 3 and Cd 4 ., re- 
spectively. Upon further cooling the system undergoes a 
first order spin-Peierls-like phase transition from a cubic 
paramagnet to a tetragonal Neel state at Tn = 12.5 K 
and 7.8 K for A = Zn 3 and Cd 4 -. The tetragonal lat- 
tice distortion induces an exchange anisotropy that lifts 
the frustration and allows the system to select a partic- 
ular spin configuration as its ground state. The nature 
of the phase transition and that of the ground state de- 
pend on the delicate balancing act between the lattice en- 
ergy cost for the distortion and the magnetic energy gain 
due to the spin ordering. Previous studies have shown 
that the ionic size of the A ion has a crucial role in the 
selection process. In the case of ZnCr204 with smaller 
Zn 2+ ions, the phase transition involves a tetragonal con- 
traction along the c-axis with I4m2 symmetry 5 - and com- 
mensurate Neel ordering 3 -. In the case of CdCr204 with 
larger Cd 2+ ions, however, the transition yields a tetrag- 
onal elongation along the c-axis with I4i/amd symme- 
try and an incommensurate Neel ordering 4 -. Their mi- 
croscopic mechanisms are not to be understood yet. A 
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FIG. 1: (Color online) Phonon spectra measured at 10 K 
from a single crystal of CdC^CU along Q=(6,8,4)+(h,-h,l) 
direction. Each solid (red) line shown is a fit to Gaussians. 



theory based on Dzyaloshinskii-Moriya interactions was 
proposed to explain the static spin-lattice coupling in 
CdCr 2 04^. This theory, however, does not provide an ac- 
curate account of the one-to-one correspondence between 
the tetragonal distortion and the Neel state that were ex- 
perimentally observed, and the microscopic mechanism of 
the static spin-lattice coupling is yet to be understood. 

Another issue is whether or not the spin-lattice cou- 
pling is also dynamic in nature. Recently, reflectiv- 
ity data obtained from samples of ZnC^O^ 7 -"— and 
CdC^O^^r— have been presented as an evidence for the 
spin-phonon coupling in these materials. For CdCr2 04, 
one of the modes around 45 meV splits into two modes 
while the other three modes do not, below Tn- The 
45 meV mode is a F\ u symmetry mode that domi- 
nantly involves displacements of the magnetic Cr 3+ ions. 
Throughout comparing the intensities of the two split 
modes, they concluded that the higher energy mode of 
the two is doublet and the lower energy mode is sin- 
glet. The c-axis elongation in the tetragonal phase of 
CdCr204, however, should yield a lower energy doublet 
mode (yz and zx) and a higher energy singlet mode (xy) . 
This led them to conclude that the phonon anomalies 
could not be simply explained by the tetragonal dis- 
tortion rather than spin-phonon coupling must be in- 
volved 1 ^. However, the reflectivity spectra is a limited 
technique in studying the lattice dynamics because it is 
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FIG. 2: (Color online) Phonon dispersion relations as a func- 
tion of Q and Tiu) obtained at 10 K CdCr2 04. Red diamonds 
indicate the experimental values while black lines show cal- 
culated dispersions, (a) Along Q=(4, 4, h), (b) along Q=(h, 
0, 0), (c) along Q=(6, 4-h, 4+h) and (d) along Q={8, -h, h). 
Here Q = f + q where f is a reciprocal lattice vector. 



a probe of phonons at a zone center (Q = 0). 

Here we report our inelastic synchrotron x-ray scatter- 
ing (IXS) measurements on a single crystal of CdCr204. 
We characterized the lattice vibrations of CdCi'204 both 
in the cubic Fd2>m and the tetragonal IA\/amd phases 
by mapping out the phonon dispersion relations up to 50 
meV along high symmetry directions. Our analysis of the 
data done using a computer software, GULP—, leads us 
to conclude that the energy shifting of the phonon modes 
around a zone boundary is simply due to the ordinary ef- 
fects of the lowering of the crystal symmetry. 

A ~ 100 mg single crystal of CdCr 2 04 (space group 
Fd3m, a=8.58882 A for T=10 K) was grown using a 
flux metho d 12 ' 13 . The IXS measurements were done at 
the high resolution BL35XU beamline at the SPring-8 in 
Japan. Using the backscattering geometry allows obtain- 
ing a large angular acceptance for both the monochroma- 
tor and the analyzer and an excellent energy resolution 
of ~ 1.6 meV over the entire range of energy up to ~ 50 
meV. A Si(ll 11 11) monochromator was set to produce 
incident x-ray with Ei = 21.747 keV onto the sample. 

We have first characterized the lattice vibrations in the 
cubic phase. Fig. [JJ shows the phonon spectra measured 
at several different wave vectors along the (h,-h,l) direc- 



tion in the cubic phase. Similar measurements were done 
along the (h,0,0) and (0,-h,h) directions. A primitive cell 
of the cubic (Fd3m) crystal structure has two formula 
units of CdCr204 and so 14 atoms. Each atom generates 
three vibrational modes, thus there are 42 phonon nor- 
mal modes in total. Thus, we have fit the data with 42 
gaussians which correspond to 42 possible phonon modes 
with a fixed energy resolution. The results of the fits are 
summarized in Fig. [2] as a function of momentum and 
energy. 

The dispersion relations of the phonon modes have 
been analyzed using group theory. Theoretical group 
calculations for Fd3m crystal structure yield 39 optical 
modes at a zone center (r point). The optical modes are 
made of the following symmetries, 

r = A lg + E g + F lg + 3F 2g + 2A 2u + 2E U + 4F lu + 2F 2u , 

(1) 

where the symbols. A. E, and F, represent singlet, dou- 
blet, and triplet, respectively. Among them, the Ai g ,E g , 
and Fig modes are Raman active, the F\ u modes are in- 
frared active, and the rest are inactive in both Raman 
and infrared reflectivity measurements. The dispersions 
of the phonon modes were calculated by the computer 
software package, GULP—, with a rigid ion potential 
model. In this model, the potential energy for the lattice 
vibrations can be written as 

where K;/< is a short range force constant, (ru> — tq , ) is 
the displacement vector from the equilibrium position of 
the Z-tli and T-th atoms and Z; is an effective charge on 
the i-th atom. The first term of eq. ^ corresponds to the 
short range repulsive forces, while the second part is due 
to the long range Coulomb interactions. For simplicity, 
the Coulomb interactions are considered only between 
the nearest neighbor atoms. In the calculations we used 
six short-range force constants described in Fig. [3] (a) . 
An initial set of those constants was used to generate a 
reasonable fit to the experimental data, and the experi- 
mental data were re-fitted to 42 gaussians to obtain the 
peak positions, i.e., characteristic energies of the phonon 
modes at different Qs. Because some of the modes are 
too close in energy to be distinguished by the experimen- 
tal resolution, the number of peaks was reduced to 16 or 
20 depending on which Q when the data were refitted 
to 42 Gaussian. We repeated this process several times 
until the comparison between the calculated dispersions 
and the experimental data were converged. The final op- 
timal values of the force constants are listed in Table U 
and the resultant dispersion relations are shown by the 
black lines in Fig. [2] 

Table [Til lists the energies of the phonon modes at a 
r point. The values of energies arc consistent with the 
previous phonon studies on other spinel compound o 14 ' 15 . 
Here, at a T point the (i^u and Fi u ) modes with < 
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FIG. 3: (Color online) Six short-range force constants used 
for a rigid ion model phonon calculation described in the text 
(a) and the atom vibrations in CdCr2C>4 phonon modes at T 
point (b-f) in the cubic phase, (b) Cd vibrations dominant 
mode of F2 9 at 16.3 meV. Cr vibrations in (c) F2u at 22.1 
meV, (d) E„ at 39.1 meV, and (e) Fi« mode at 44 meV. 
(f) O vibrations in Ai 9 mode at 87.1 meV. Black arrows are 
directions of the vibrations. 
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FIG. 4: (Color online) The IXS phonon spectra as a function 
of energy, hu;, obtained from a single crystal CdCr2 04 along 
<2=(6,8,4) + (h,-h,l) at 10 K (black diamond) and 6 K (red 
asterisk), in the energy range (a) from 2 to 50 meV, (b) from 
2 to 20 meV: Cd dominant modes, and (c) from 20 to 50 meV 
: Cr dominant modes. 



TABLE I: The values of the short-range force constants and 
effective dynamical ionic charges 

Force constants Values Effective charges Values 
(eVA) (e) 



Cd-O (Ki) 


11.22 


Cd 


0.44 


Cr-0 (K 2 ) 


9.77 


Cr 


1.64 


Cr-Cr (K 3 ) 


0.98 


O 


-0.93 


0-0 (K 4 ) 


0.46 






0-0 (K s ) 


0.70 






0-0 (K 6 ) 


0.00 







TABLE II: Calculated phonon energies at a zone center above 
and below Tn in CdCr2 04. When the system transits from 
cubic to tetragonal, most phonon modes change only by a few 
tenth percent of their energies, consistent with the few tenth 
percent change in the atomic distances (hui cx r _3//2 ). The 
oxygen-stretching Ai 9 and A2 U modes, however, change by ~ 
2 % in energy. This can be understood by the fact that the 
atomic distances between oxygen atoms reduce by 0.8 % (e.g. 
2.596 A in cubic to 2.576 A in tetragonal phase) so that the 
energy changes by 1.3 % 



Fd3m, 10 K 


lii/amd, 6 K 


Species 


E (meV) 


E (meV) 


Species 


Ai 9 


87.1 


85.3 


Ai 9 


A 2u 


60.0 


60.0 


Bi„ 


89.9 


88.1 


E s 


49.4 


48.2, 49.1 


Alg + Big 


E u 


39.1 


39.1, 39.3 


Ai„ + Bi u 


66.4 


65.6, 66.3 


Fi 9 


54.0 


54.1, 54.2 


A2 9 + E 9 


F 29 


16.3 


16.2, 16.4 


B 29 + E 9 




63.5 


63.5, 63.7 






77.5 


76.0, 77.5 




Fi„ (TO) 


21.6 


21.1,21.2 


A 2u + E„ (TO) 




44.4 


44.3, 44.4 






58.6 


58.3, 58.7 






75.2 


74.2, 75.1 




Fi„ (LO) 


21.7 


21.7 


A 2u + E u (LO) 




44.8 


44.6 






69.6 


68.6 






85.3 


85.3 




F 2u 


22.1 


22.3, 22.3 


B2u + E„ 




46.0 


46.0, 46.1 





fiui < 22 meV involve mainly vibrations of the heavy Cd 
ions, the (F lu , F 2g , and E u ) modes with 22 < Tllo < 
45 meV involve mainly vibrations of the Cr ions, while 
the modes with hu > 45 meV involve mainly O ions 
except modes with Tiuj ~ 60 meV which have dominant 
Cr vibrations. Thus, the (F\ u , i^g, and E u ) modes are 
expected to show anomalous behaviors when the system 
undergoes the phase transition if spin-phonon coupling is 
involved. Eigenstates of some typical vibrational modes 
at the zone center are illustrated in Fig. [3] 

After the careful characterization of the lattice vi- 
brations in the cubic phase, we have also studied the 
tetragonal phase. First, we confirmed the symmetry 



4 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 



h h 

FIG. 5: (Color online) The peak positions of phonon modes 
(circles) and calculated dispersions as a function of g=(h,-h,l) 
at 10 K (black solid lines) and 6 K (red dashed lines). The 
size of each circle is proportional to the integrated intensity 
of each peak, (a) and (b) show the dominantly Cd and Cr ion 
vibrational modes respectively. 

of the tetragonal structure to be IAi/amd not 74i22 
by finding out that elastic nuclear reflections are ab- 
sent at (6,0, 0) cubic = (3, 3,0)totra and (8,2,0) cubi c = 
(3, 5, 0)tctra- Second, we have performed IXS measure- 
ments to study lattice vibrations in the tetragonal phase. 
Fig. [4] (a) shows the data obtained upto 50 meV at 6 K 
(< Tn), over-plotted with the 10 K data of the cubic 
phase at three different wave vector positions along q = 
(h, -h, 1). As shown in Fig. 2] (b) and [5] (a), among 
the mainly Cd vibrating modes upto 20 meV, two modes 
show changes between the two temperature data: for Q 
= (6.1, 7.9, 5) the peak centered around 18 meV becomes 
broad at 6 K and for Q = (6.38, 7.62, 5) the peak cen- 
tered at 12.0 meV at 10 K shifts to 11.4 meV at 6 K. In 
Fig. 2] (c) and|S](b) which show the mainly Cr vibrating 
modes, the peak that appears at 32.1 meV at 10 K moves 
to 31.0 meV at 6 K for Q = (6.1, 7.9, 5). At the same 
Q position, the two peaks that appear at 45.8 and 47.2 
meV at 10 K move to 44.8 and 46.3 meV below T N . The 
other peaks do not show any obvious difference between 
two temperatures. Similar measurements were also per- 
formed along Q = (7,7,5) + (h,h, 0.2), and only those four 
modes seem to shift down by ~1 meV when the tempera- 
ture changes from 10 K to 6 K. Now, the question is that 
the observed changes in phonon spectra are simply either 
due to the cubic-to-tetragonal structural distortion or a 
spin-phonon coupling. In order to address this issue, we 
used the rigid ion model to reproduce the phonon spectra 
in the tetragonal phase. 

There are 39 optical modes in the IAi/amd tetragonal 
phase at a zone center, which can be described as 

T =2A lg + 2A lu + A 2g + 4A 2u 

+ B lg + 4B lu + 3B 2g + 2B 2u + AE g + 6E U . (3) 



Therefore, the triply-degenerate modes split into a sin- 
glet and a doublet mode in the tetragonal phase, and the 
doubly-degenerate modes split into two singlet modes. 
In the cubic phase, the distance between nearest neigh- 
bor Cr ions is r^N = 3.037A. In the tetragonal phase, 
fNN = 3.034A is for the bonds in the afr-plane and 
i~nn — 3.040A is for the bonds in the etc- and 6c-planes. 
Since the difference in r^N is small, we assumed that the 
force constants, K;p, remain unchanged after the tetrag- 
onal distortion. The effect of the lattice distortion is 
reflected as the changes in Coulomb forces. The calcu- 
lated energies of the phonon modes at the T point for 
the tetragonal phase are listed in Table [TT1 and their 
dispersions are shown by dashed (red) lines in Fig. [5] 
Even though the simple rigid ion model does not per- 
fectly reproduce our data, the calculation shows that 
the four modes around 13.2, 14.8, 33 and 48 meV at 
q = (0.1, —0.1, 1) position in cubic phase shift to 12.8, 14, 
32.5 and 47 meV in tetragonal phase. This is consistent 
with our experimental observation at the same q point 
where the modes with Klo ~ 13, 32, and 46 meV at 10 K 
shift to lower energies by ~ 1 meV at 6 K. The observed 
changes in our x-ray phonon spectra seem to be due to 
the ordinary mode shifting by the lowering of the crystal 
symmetry, rather than the spin-phonon coupling, which 
was previously claimed by reflectivity measurements^— . 
It remains to be seen whether or not a zone-center and 
a zone-boundary should behave differently in this com- 
pound. 



In conclusion, we have searched possible phonon 
anomalies that might be associated with the origin of the 
phase transitions in the frustrated magnet CdCr 2 04 us- 
ing inelastic synchrotron x-ray scattering. This was to 
address the issue of whether or not a static as well as 
dynamic spin-lattice coupling is part of the mechanism 
leading to the phase transition. Our results show that a 
phonon mode energy around a zone boundary with Huj ~ 
45 meV shifts down as the crystal structure changes from 
cubic to tetragonal. The shifting can be well reproduced 
by a rigid ion model with the same force constants as 
used for the cubic phase. 
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